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In this paper, we discuss and mathematically compute the eigenvalues and 
the characteristic polynomials of special even square matrices of orders 
4x4 and 8x8. Also, we introduce two 8th order compound magic squares. 
The computed values are verified using Maple software. First, for the 
4th order square matrix, the characteristic polynomial was derived to be: 
A(A — 2s)(A? + 40) with the eigenvalues: 0,2 s, and two other conjugates. In 
further analysis, we performed numerical classification of the squares for the 
matrices of order 4. Second, for the 8th order magic square, the 
characteristic polynomial was obtained in the form: 43(A — 4s)(A* + 0A? + 0) 
where Q, 0 are constants; the eigenvalues are 0,4 s, F/M, Fy à2; where Aj, Az 
are the roots of the quadratic equation: A27+0A+0=0. Third, for the 
franklin square, we obtained the eigenvalues 0,4 s, and the roots of the 
equation: 4? + ad +b. Finally, we suggested a hybrid image encryption 
technique based on Franklin magic square matrices and improved 
substitution technique. The proposed a grayscale image encryption/ 
decryption algorithm uses circular rotation of bits and Franklin magic 
squares’ properties in conjunction with substitution techniques to obtain a 
very secure algorithm against attacks. 
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1. INTRODUCTION 

In this investigation, we derive the characteristics polynomial and compute the eigenvalues related 
to special matrices called magic squares. Also, a semi magic square is a matrix that has the same sum of the 
matrix entries of every row and column, called the magic constant. The semi magic square has the sum of the 
main diagonals equal to the magic constant. On the other hand, the n™ order natural magic square is a square 
matrix with integer entries from one to n?. The magic constant is obtained as (1). 


=n(n? +1) (1) 


The combinations which appear in the columns, rows, and both diagonals of this square are the only 
distinct three elements’? mixture of the numbers from 1 to 9 that sum up to the value of 15. A self 
complementary magic matrix of the order 2k+1 that has the (2k+1)s as a magic constant is still considered a 
magic square satisfying the conditions (2). 


QAlk+1)(k+1) = S & Qij + Agok+2-)2k+2-j) Z 2S,1 Si <Sk,1 Sj <2k+1 (2) 
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Also, a 2k self complementary magic matrix with 2ks as a magic constant is a magic square having 
the condition (3). 


aij + Q(2k+1-i)(2k+1-/) = 2s,1 <i< k,1 <j < 2k (3) 


In other words, a natural self-complementary magic square is a magic square of order n such that the sum of 
both elements of each couple of dual (opposite entries) is usually equal to n° +1. 

An off-diagonal is a combination of two parallel diagonal lines to the same main diagonal. The two 
parallel diagonal lines must occur on opposite sides of the primary diagonal and they can only be combined if 
the combination has the same number of entries as the primary diagonal. Two examples of an off-diagonal 
line are 10, 3, 8, 13 and 7, 14, 9, 4 are shown in Figure 1. The Figure shows 3 off-diagonals corresponding to 
the main diagonal. The square in Figure 1 is called semi pandiagonal since the off-diagonal lines 10, 3, 8, 13 
and 7, 14, 9, 4 sum to the magic constant. A pandiagonal square is a magic square where the off-diagonals 
sum to the magic constant. It must be noted that even by swapping all rows and columns of a pandiagonal 
square, the matrix remains pandiagonal. 

The results of different order squares are verified using linear algenbra while Maple software is used 
to pefrom lengthy calculations; some complicated calculations are difficult to manipulate manually and must 
be performed using a computing device. Additionally, we introduce special cases of 4 by 4 squares which 
have a simplified form of the eigenvalues. Also, we give some numerical examples that substantiate our 
results. 


1 7 10 16 
4122 5 3 
s$ 2 I 9 
ll 13 4 #6 


Figure 1. Natural semi pandiagonal square 


2. LITERATURE REVIEW 

There are many papers in literature related to magic squares and their applications. Nordgren [1] 
proved using Mattingly well-known theorem that magic squares called quasi-regular (QR) have signed 
eigenvalues pairs are similar to that of regular magic squares. The researchers concluded that the odd powers 
of the QR magic matrices are magic squares that can be obtained from the QR condition. They generalized 
that all 4" order pandiagonal magic squares are QR because they are all most-perfect (MP). Also, the authors 
argued that the QR conditions holds for all 5" order pandiagonal magic squares but the same cannot be said 
of higher order matrices. Similarly, Stephens [2] studied 4" and 5“ order magic squares and elaborated on the 
previously obtained eigenvalues. An n by n matrix A = (A; ;) is called a pentadiagonal matrixif A; ; = 0 
when |i — j| > 2 condition is true. Mattingly [3] argued that even order regular magic squares are singular 
because they always have a zero eigenvalue and since each positive eigenvalue has the same algebraic 
multiplicity as the negative one. Chu et al. [4] analyzed magic matrices for which the moore-penrose inverse 
is also magic. They introduced singular magic matrices in which the numbers in the rows and columns and in 
the two main diagonals all add up to the same sum. They then introduced the concept of a “philatelic magic 
square” as a square arrangement of images of postage stamps for which the nominal values make a magic 
square. Ibrahim and Salman [5] presented definitions and some special properties of magic squares with 
proofs of two theorems. Neeradha and Mallayya [6] presented the famous strongly magic square, which has 
the strong property that the sum of the entries of the sub-squares taken without any gaps between the rows or 
columns is also a magic. The authors detailed the properties of 4x4 strongly magic squares dot products and 
different eigen values and eigen vectors properties. Zhaolin et al. [7] suggested a new Sylvester-Kac matrix, 
i.e., Fibonacci-Sylvester-Kac matrix where they discussed the eigenvalues, eigenvectors and characteristic 
polynomial of the matrix in two categories based on whether the Fibonacci-Sylvester-Kac matrix order is odd 
or even. Additionally, the authors gave the explicit formulas for its determinant and inverse. Abdulkhaev and 
Shirokov [8] discussed the characteristic polynomials in (Clifford) geometric algebras Gp, q of vector space 
of dimension n=p+q. They presented explicit formulas for all characteristic polynomial coefficients in the 
case n=5. The formulas involved the operations of geometric product, summation, and conjugation. 

Tridiagonal and pentadiagonal matrices appear in several areas of mathematics and engineering; 
specially the ones that involve linear systems of differential equations. Since every 3 by 3 matrix is a 
pentadiagonal one, it is clear that not every pentadiagonal matrix is comparable to a symmetric matrix. 
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Alvarez et al. [9] considered two classes of pentadiagonal matrices and obtained recursive formulas for the 
characteristic polynomials and explicit formulas for the eigenvalues of these classes of pentadiagonal 
matrices. They showed that if A is a pentadiagonal matrix and if it satisfies some certain conditions on the 
sign of the product of the entries, then A is similar to a symmetric pentadiagonal matrix. Thus, all the eigen 
values of A are real and A is diagonalizable. Elouafi [10] explained that the characteristic polynomial for 
such matrices is the product of two polynomials presented in terms of Chebyshev polynomials. 

Al-ashhab [11] considered some properties of the magic 4 by 4 squares such as their determinant; 
while Al-ashhab [12] introduced new special types such as the four corner magic squares. The researcher 
explained ways of counting them using parallel programming. He also discussed the characteristic 
polynomial of Franklin squares. Rosser and Walker [13] described the general structure of pandiagonal 
magic 4 by 4 squares. The authors proved for the first time a well-known structure. Ahmed [14] described 
many features of Franklin squares, which is used in this paper in a cryptographic application of magic 
squares. Fahimi et al. [15] classified twelve groups of 4 by 4 squares where he discussed these groups with 
the minimum values of electrostatic potential. The potential and its related properties are calculated on the 
grids numerically using MATLAB 2018a. Equipotential points and certain constants are found among the 
electrostatic potential sums along horizontal and vertical lines on the square lattice. 

Rungratgasame et al. [16] introduced new special magic squares and called them reflective magic 
squares, corner magic squares, and skew-regular magic squares by combining the concepts of magic squares 
and linear algebra. They found the dimensions of the vector spaces of these magic squares under the standard 
addition and scalar multiplication of matrices using the rank-nullity theorem. 

Many algorithms and schemes have been proposed to protect data of various formats including 
images. Ozturk and Sogukpinaar [17] classified image encryption techniques into three types: position 
permutation, value transformation and visual transformation. For example, Caesar's magic square works by 
removing spaces between words of sentences. The letters are arranged in a square starting from the top left 
downward, then moving to the second column from the top and so on. It is preferable to have the letters of 
the sentence as squares of an integer, such as 1, 4, 9, 16, 25, 36, and so on, so that the square root of the 
letters can be obtained. Al-qdah [18] proposed and implemented a secure two operation image encryption 
technique by dividing a selected image into some blocks before applying bit rotation using one random key 
followed by chaotic map scrambling using another random key. Gupta et al. [19] introduced a hybrid image 
encryption technique for protection against statistical attacks using Arnold map and S-box of advanced 
encryption standard (AES). The Arnold cat map spatially shuffles the pixels of the image and then the image 
is encrypted using substitution-box (S-box). Irom and Ningthoujam [20] considered the applicability of a 
magic squares/weak magic squares matrix of any order in evaluating numerals for encryption and decryption. 
The authors indicated that singly even n weak magic squares can produce possibly different ciphertext from 
plaintext than that of the actual magic squares. Also, they added dummy letters to the English letters to 
eliminate the duplication of vowels in encrypting a message. 

Alattar and Rahma [21] proposed a hybrid method to encrypt images using franklin magic square 
matrices and improved substitution techniques. They developed two image or text encryption algorithms 
based on 5 by 5 Magic Squares using GF(P) and GF(2°): the first used a message length = 10 and the second 
used a message length = 14. A number of rounds were added and a mask is used in the even round for an 
addition operation while the odd round is used for multiplication operation so that the text resulting from the 
first round is used as input text for the next Round. Al-Hashemy and Mehdi [22] presented a hybrid image 
encryption algorithm based on chaotic systems and image squares. The encryption keys are generated 
chaotically equal to the size of an image, which is then subdivided into equal sub images; each sub-image is 
multiplied by a magic matrix. The resulting matrices are XORed produce one encrypted image. 

Alattara and Rahma [23] developed an encryption algorithm using magic squares (5 by 5) with multi 
message lengths for added security. The algorithm used two rounds and selected a message interchangeably. 
The key was fixed in one position while the remaining positions are filled with the message; some sums were 
calculated to obtain an encrypted text. The authors concluded that their proposed algorithm gave better 
encryption speeds and added levels of security. Umar [24] presented an added security layer to the Rivest- 
Shamir-Adlemen (RSA) algorithm approach using a magic square of order 32 to generate a non-duplicate 
random numbers to represent the numeral aspect of the message instead of the ASCII values and also the 
author used the magic squares values as a key in the encryption process. Mohammed and Hasan [25] 
suggested an encryption technique that hides a message in a submatrix of order 4 by 4, which is selected from 
a 16 by 16 magic square matrix. The author used two stages of hiding the ciphertext using a magic square of 
size 3x3 and Latin square of size 3x3 for more permutations and to obtain an inverse matrix for the 
decryption operation. The algorithm hides the ciphertext into a 16 by 16 matrxi with 16 sub matrices so that 
duplicates are eliminated. The elelments of the matrices were polynomial numbers of a finite field of degree 
Galois fields GF (28). 
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3. METHODS AND DISCUSSIONS 

In this section, we particularly discuss the charactersitics of magic squares of sizes: 4 by 4, 8 by 8, 
nested 6 by 6, and franklin squares. We first perform some mathematical proofs and calculations before we 
utilize Franklin magic squares to develop a cryptographic technique in subsection 3.6. 


3.1. Magic square (4 by 4) 
It is accepted that the magic constant is always an eigen value since one of the corresponding eigen 
vectors has all ones. Hence, any 4 by 4 semi magic square has its characteristic polynomial in the form (4). 


(A — magic constant)( A? +a 2? +ba +c) (4) 


The determinant of the pandiagonal magic square 4 by 4 is zero. Generally, the pandiagonal magic square of 
size 4 by 4 has the structure shown in Figure 2. 


A B C 2s-B-C-A 
E 2s-B-A-E A+E-C B+C-E 
=C A+B+C-s s—-A s-B 
C+-A-E +C-B+E s-E A+B+E-= 


Figure 2. Symbolic pandiagonal square 


The characteristic polynomial of the magic square was calculated using Maple software to be (5): 
A(A — 2s)( 4? + 40) (5) 


where 


@ = C? — A? — s? — AE —-2BE—CE-—AB+BC+s(2E+2A+B-C) (6) 


we can see that the characteristic polynomial has two missing coefficients (a and c) and can be rewritten as 


(7): 
(A — 2s)(A3 + 404) (7) 


Mathematically, we can obtain the characteristic polynomial in (1) using the trace eigen value theorem which 
postulates that the sum of the eigen values is the trace (in this case 2s). We can also conclude that both eigen 
values must be conjugates since the magic square has 2s as eigen values, and it has been proven that its 
determinant is zero (0 is the second eigen value). Figure 3 shows another example of a 4 by 4 magic square. 
The figure presents a semi magic square of order 4 that has the characteristic polynomial with no missing 
terms (8). 


(A — 60)( 4? — 120 4 ? — 3453 A + 980) (8) 


Figure 3. A 4th order semi magic square matrix 


Therefore, we can generally conclude that semi magic squares don't have have any missing 
coefficients in their characteristic polynomial. On the other hand, the standard magic square in Figure 4 has 
2s as a magic constant. We observe that the matrix has seven independent variables in addition to the s 
variable. The characteristic polynomial, which was calculated by Maple, has one missing coefficient (a=0) 
given in the form (9). 


(A — 2s)(43 +bA+c) (9) 
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p j 2s-j-h-p 
v q 2s-q-f Êv 
6s-2h-j-29-Fv-p 2sq fħy2s 2hH+lq-4s+v+p 
f+htj+2q-4s 25-j 4s-h-2q-f 


Figure 4. Magic square matrix with constant 2s 
The mathematical justification for the above form can be performed using the trace coefficient 
theorem which states that the trace (in this instance 2s) is the negative of the second coefficients; then the 
characteristic polynomial must have the form (9). 
(A — 2s) (23 +b A+ c)=A* — 2s2? + bd? + (c — 2bs)A — 2cs 


(10) 
In general, semi pandiagonal magic squares with magic constant 2s has six independent variables in addition 
to the s. Therefore, the square matrix can be written in the form given in Figure 5. 


f d z—d+t+m-f 
z t+m-f 2s-2t-m+f-z+d-g eit 
m g z-d+t 2s—2-t+-m—ztd 
t 2s—t-m-¢ zH+m-fd 


Figure 5. Semi pandiagonal magic square 


Furthermore, we calculated the semi pandiagonal magic square using Maple software; it has the 
characteristic polynomial in the form: 


(A —2s)(2 +b A+ c0) 


(15) 
Figure 6 shows an example of a semi pandiagonal magic square with two independent variables. The 
characteristic polynomial was calculated as (16). 


(A — 2s)(8—(14f + 10s — 2fs + f? — 15) A — 2M) 


(16) 
Where: 
M = 4s? + (2f — 34)s? + (79 — f? — 14f)s + 2f? + 26f — 36 (17) 
23-3 f + -6&1 
2 l-f 2:4 l 
1 5 -2 25-4 
0 23-6 2-f f+4 


Figure 6. Special semi pandiagonal magic square 


in some particular cases, a semi pandiagonal magic square with less free elements can have a simple form of 


the eigenvalues; for example, Figure 7 shows a special case with five free elements. It can be observed that 
the eigen values are independent of d and are actually in the following form (18). 


0,2s,+2,/(s —t—z)(m+t-—s) 


(18) 
25-m-t-= 5-d d mHz 
z dhm=H 2st sz 
m d+s-t-z td+z zm 
t s—damn+z dHn-z 


a-t} 


Figure 7. Special semi pandiagonal magic square type 1 
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Figure 8 illustrates one more special case where the last two eigen values are indpenedent of s and j. 


a-bjts b sb b-t 


J 5-£ E =—] 
2b ets -g  sj-łb 
=b- -b bH ] 


Figure 8. Special semi pandiagonal magic square type 2 


Specifically, the eigen values are: 0,2s, +24/ b? — g?, Figure 9 shows an additional special case where all the 
eigen values can be real numbers. 


Zs f d —Lf 
0 -Ë dels zd 
0 z -d d-ztis 
O 3s E-fd d+f 


Figure 9. Special semi pandiagonal magic square type 3 


The eigen values are given in the form (19). 

2s,2s,-st/f?+2fg —3g9’?+ 4dg + (6g — 2f — 4d)s — 3s? (19) 
The discriminant of the quadratic quantity under the root sign with respect to s is (20). 

16(df + Pœ + f?) > 0,Vd,f #0 (20) 


Hence, if f and d are both nonzero then there is always a bounded interval of s (between the roots of the 
quadratic quantity) that has four real eigen values regardless of the g value. 


3.2. Magic square 8 by 8 

We construct an 8 by 8 pandiagonal magic square using four separate pandiagonal magic squares of 
sizes 4 by 4 with each having the same magic constant. The square matrix structure is displayed in Figure 10. 
We require for the four submatrices to become pandiagonal that some conditions on D, F, H, L, N, O and P 
and the corresponding small letters hold. As a result, the generated 8 by 8 square becomes of magic type; but 
in order to obtain a pandiagonal 8 by 8 square as shown in Figure 11, we must impose more conditions on the 
independent variables: At+i=a+I, B+j=b+J, C+k=c+K, E+m=e+M. This correspondence will be passed onto 
the dependent variables automatically, e. g. D+l=d+L. 


A B c D I J K L 
E F G H M N (0) P 
s-C -D :-A s-B s-K -L s-I -J 
-G :-H -E -F s-O s-P s-M s-N 
a b c i j k 1 
2 f gZ h m n o p 
s-c s—d s-a s-b s-k -l1 5-1 53- 
5-g s-h s-e s-f 5-0 s-p s-m s—n 


Figure 10. Special pandiagonal magic square of size 8 by 8 


The matrix rank is five and its characteristic polynomial is (21). 


X3 (A — 16)(A* — 38402 + 55328) (21) 
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Figure 11. Special numerical pandiagonal magic square of size 8 by 8 


Thus, we conclude that its rank is generally five and expects that the characteristic polynomial to have the 
form (22): 


BA — 4s)(A* + OA? +O) (22) 


where Q, @ are constants. Also, when we perform thesubstitution x = A”, we obtain the roots of this equation 
(the eigenvalues): 0, 4s, F/M, Fy Az, where A, à; are the roots of the quadratic in (23). 


x?+0x+0=0 (23) 
The proof of this eigen values conjecture is achieved by computing the determinants using a 
computing device since the amount of computation of the determinants is extremely huge. Now, we again 


consider the characteristic polynomial and a special type of matrices that have the form (26). 


WA — 4s)? +o) (24) 


Where ¢ is a constant. We arrange two 4 by 4 pandaigonal magic squares A and B (magic sum 2s) that satisfy 
the conditions in the following matrix (25). 


v= (G) os 


Accordingly, we conclude that the nullspace of this 8 by 8 matrix contains the following four independent 
vectors (x;) such that (26). 


xi = —1, xiş4 = 1, x; = 0 forj + i i + 4,i = 1,2,3.4 (26) 
Thus, the characteristic polynomial must be in the form (27). 
24 (A — 4s) (2? + bà + c) (27) 


But the square of the roots of this cubic polynomial must be eigen values for D?. The structure of D? is as 
(28): 


F a) (28) 


where G and H are matrices, which has the following structure (29). 


P ọ of 

|X T DE 
N=\a $ po (29) 

V X T 


In addition, we observe that G and H have the same nullspace. If we take the vector (v1, V2,, V3, v4)! in this 
common nullspace, then a basis of the nullspace of D? can be chosen as (30): 


(V4, V2,, V3, V4, 0,0,0,0)*, (xi), xi = —1, X44 = 1, xj =0 forj #i,i+ 4,i = 1,2,3.4 (30) 


We notice that N multiplied by the vector (1,1,1,1)' yields (8s?,8s?,8s7,8s)*. On the other hand, N 
multiplied by the vector (1,—,1,1,—1)* yields (t, —t, t, —t)t for a real number t. Therefore, D?is a magic 
square with a magic sum 16s”; and this magic sumis one of the Deigen values. The eigen vectors are 
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multiples of the vector one. Also, due to these properties we can deuce that the value t, + t is an eigenvalue 
of D? where: 


tı = G(1,-1,1,-1),t. = H(1,-,1,1, -1)', (31) 


an eigen vector can be selected as (ty, —t,, ty, —t,, tz, -t>, t2, -t,)*. Actually, we can check that there exists 
a second linear independent eigen vector for this eigen value, and 


X(X — 16s?)(X — (tı + t2)) (32) 


is the minimum polynomial of D?. Therefore, the only eigen values of D? are 0, 16s? and ty + tz. The other 
eigen values of D (except 0 and 4s) must sum up to zero and their squares is one of the three eigen values of 
D?. Thus, the eigen values must be 0, 4s and Fyt; +t, and the characteristic polynomial is given by (33). 


A5(A — 48)(A? — (ty + t2)) (33) 


3.3. Franklin squares 

Franklin squares are semi magic squares of order 8. They have some similarities to pandiagonal 
magic squaresof size 8 by 8 and resemblance to the special type of 8 by 8 square, which is a combination of 
four squares. As explained by Franklin, each row and column of the square has the common sum 2s and half 
of each row or column sums up to half of s. In addition, each of the "bent rows" (as Franklin called them) 
have the sum 2s. Franklin described the 8 by 8 magic square and its special properties. Figure 12 shows an 
example of a Franklin magic square with magic sum 260. 


52 61 4 13 20 29 36 45 


14 3 62 51 46 35 30 19 
53 60 5 12 21 28 37 Gad 
11 6 59 34 43 38 27 22 
55 58 7 10 23 26 39 42 
9 8 57 56 41 40 25 24 
50 63 2 15 18 31 34 47 
16 1 64 49 48 33 32 17 


Figure 12. Numerical Franklin magic square 


The general form of a Franklin square with magic constant 2m is (34): 


D= (i (34) 


where A, B, C and D are semi magic 4 by 4 squares. Neeradha and Mallayya [6], it was proven that the 
characteristic polynomial is in the form (35): 


25 (A — 4s)(A2 + ad + b) (35) 


where the constant a depends only on the value of two entries in the square. Actually, there are only eight 
independent variables in addition to the magic constant. 


3.4. Compound magic squares 8x8 

In this section we present an 8" order square shown in Figure 13. The matrix consists of four 
nonregular semi magic squares where each one has the number 130 as a magic constant value. We are going 
to show that the characteristic polynomial is simpler due to the form and the way its constructed. 


1 62 44 23 16 1 26 37 
63 4 22 41 30 13 40 27 
24 43 61 2 25 38 15 2 
42 21 3 64 39 28 49 14 
17 46 60 7 32 35 10 3 
47 20 6 57 34 29 6 11 
8 59 45 18 9 4 31 36 

8 19 48 5 12 33 30 


Figure 13. Compund magic square that can be divided into four 4th order semi magic squares 


Indonesian J Elec Eng & Comp Sci, Vol. 27, No. 3, September 2022: 1479-1489 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1487 


The matrix in Figure 13 has a nontrivial nullspace and one of its bases is given in (36). 


(1,0,1,0, 1, 0, 1, 0), ( 1, 1,0,0,1,1,0,0), ( 1, 0, 0, =1, —1, 0, 0, 1 ). (36) 
Also, the characteristic polynomial is given as (37). 
(X — 260)X*(X? + 8X? — 5120X + 81920) (37) 


Therefore, the nonzero eigenvalues are 260, 17.53, 56.77 and —82.3. But, there are various examples of 
squares that consist of four self complementary semi pandiagonal magic squares with the nullity of exactly 
three. In some cases they are diagonalizable. 


3.5. Nested magic square 6 by 6 

The nested magic squares of order n by n can be generally written as shown in Figure 14. Where we 
set some conditions on F, G, H, N, U and V to ensure that the middle 4 by 4 square is the exact pandiagonal 
matrix that was given in figure 2. We can mathematically prove that the characteristic polynomial has the 
form (38). 


I J K L M N 
(6) A B C D -0 
P E F G H s-P 
R s-C s-D s-A s-B s- 
T -G s-H s-E s-F -T 
U s-J s-K s-L -M \ 


Figure 14. Nested magic square of size 6 by 6 


MA — 3s)(At + H23 + LA? + IÀ + K) (38) 


In general, the coefficients H, L, I and K are nonzero. On the other hand, in some cases such as the magic 
matrix shown in Figure 15, we obtain for the coefficients in the characteristic polynomial (39). 


H=0andI] = 4+ 8g- 6a- 4j (39) 


Consequently, the odd powers disappear by imposing the condition as (40). 


j=1+2g -15a (40) 
1 2 3 4 0 2 
1 E h j 8-h-j-g 3 
2 a 8-g-h-a atg-j h-atj 2 
4 4) etht-4 4-2 +h 0 
2 jet a-h-j+4 4-a atg+h-4 2 
2 2 l 0 l 3 


Figure 15. Nested magic square of order 6 by 6 


3.6. Image encryption based on Franklin magic squares and modified substitution techniques 

Many researchers have proposed hundreds of numeric and alphabetic ciphers that transform letters 
into numbers, or symbols, or transform alphabets into other forms. Recently, many intelligent algorithms 
have been suggested with many permutations and lengthy keys of millions of kilobytes. Therefore, the 
literature is full of various algorithms with varying complexities and speed. Al-qdah [18] explained that 
cryptography is a secure technique that achieves confidentiality of messages. The author elaborated that the 
term had a specific meaning for the Greek: secret organizations is provided through cryptography at a high 
level, making sure that the exchanged information is sent in a secure way that only the authorized receiver 
can access it. Nowadys, cryptography is globally used by many organizations and individuals to protect 
different kinds of stored and network data. In this section, we propose a hybrid technique to encrypt/decrypt 
grayscale images using substitution techniques and the 8 by 8 Franklin magic square matrices of distinct 
values. Magic square encryption offers the possibility of checking the encrypted image for tamperment since 
the magic condition has to hold after applying the encryption algorithm. Thus, attacking the encrypted image 
might render the matrices not magic, which would indicate some attack has been inflicted on the encrypted 
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image. Therefore, this work designs an algorithm that checks for the franklin square condition on a 
partitioned image before applying rotational and substitution techniques. The proposed encryption and 
decryption procedures are explained below. 

Proposed encryption procedure: i) generate four magic squares of sizes 8 by 8 and number them m1, 
m2, m3, m4; ii) store only the values of the 9 independant variables from each matrix (total=36) that are 
needed to regenerate the four matrices (m1, m2, m3, m4); iii) read a grayscale image of size 256 pixels (pixel 
numerical range is: 0 -255); iv) divide the image into 4 partitions (8 by 8 pixels each); v) generate a random 
number to be utilized as a key; vi) apply circular rotation of bits on each pixel in the 4 matrices based on the 
encryption key (right rotation when an even number is encountered in the key and left rotation when an odd 
number is found in the key); vii) store the image’s four rotated partitions into 4 matrices (8 by 8 each) and 
number them p1, p2, p3, and p4. The partitioned matrices are shown in Figure 16; viii) substitute the values 
of every magic square matrix for its corresponding image partitioned rotated matrix (m1 for p1, m2 for p2, 
m3 for p3, and m4 for p4); then check each substituted matrix for franklin magic condition; ix) combine the 
substituted matrices into one 16x16 matrix to create one encrypted image. 

Proposed decryption procedure: i) read the encrypted image and the encryption key; ii) divide the 
encrypted image into 4 partitions (8 by 8 pixels each); iii) store the 4 partitions into 4 matrices (8 by 8 pixels 
each) and number them e1, e2, e3, and e4. The partitions are shown in Figure 17; iv) check the partitions for 
the Franklin magic condition; v) read the 4 matrices (p1, p2, p3, p4); vi) perform opposite bits’ rotations for 
each pixel using the encryption key (left rotation with even number and right rotation with odd number);vii) 
substitute the original image matrices’ values for their corresponding encrypted image partitioned matrices’ 
values (pl for el, p2 for e2, p3 for e3, and p4 for e4); viii) combine the substituted matrices into one 
combined image and compare it with the original image. 


pl p2 el 22 
P3 p4 e3 el 

Figure 16. The partitions of an Image into four Figure 17. The four encrypted partioned matrices of the 
matrices image 


4. CONCLUSION 

We have computed and analyzed the characteristic polynomial of magic square matrices of order n 
by n for the distinctive even cases of 4, 8 and 6. We specifically demonstrated that at least one linear term for 
any magic square can be factored; but also in some cases more linear terms might be factored in the squares: 
pandiagonal 4 by 4 and Franklin magic 8 by 8. Furthermore, we examined the structure of a non-linear 
remaining polynomial for several square matrices and discovered that some kinds of polynomials possess 
zero coefficients, where we proved this fact for some special types of 6" order squares. We demonstrated that 
such types of polynomials can be written as a quadratic equation with squared variable resulting in conjugate 
eigen values. In addition, we introduced examples of 4" order matrices: pandiagonal magic squares and 
special types of semi magic pandiagonal magic square. We considerd several special types of squares of 
order 4 by 4 whose eigenvalues could be written in a simpler form compared to the general case of a magic 
square. Also, we introduced a subset of an 8" order pandiagonal magic square with a characteristics 
polynomial whose quadratic factors have even powers. In the case of nested 6 by 6 magic square, we 
concluded that the characteristic polynomial has a 4" order equation in the factorization. The odd powers are 
cancelled under some conditions; the eigen values are 0,3s F/M. , FJ Az; where A,, A, are the zeros of a 
quadratic equation. Finally, we utilized Franklin square matrix properties in conjunction with improved 
substitution technique to design a secure cryptographic grayscale image algorithm. The algorithm examines a 
franklin square matrix and its charcatrsitics otherwise a temperament is detected. 
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